Quantum critical behavior in strongly interacting Rydberg gases 
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We study the appearance of correlated many-body phenomena in an ensemble of atoms driven 
resonantly into a strongly interacting Rydberg state. The ground state of the Hamiltonian describing 
the driven system exhibits a second order quantum phase transition. We derive the critical theory 
for the quantum phase transition and show that it describes the properties of the driven Rydberg 
system in the saturated regime. We find that the suppression of Rydberg excitations known as 
blockade phenomena exhibits an algebraic scaling law with a universal exponent. 
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The concept of universality is a powerful tool for the 
understanding and characterization of complex phenom- 
ena in different fields of physics. The most pronounced 
example represents the universal scaling in systems close 
to a second order phase transition and the characteri- 
zation of the transition in terms of critical theories and 
universality classes Its main assertion is that in the 
presence of a diverging length scale, physical observables 
become independent on the precise microscopic realiza- 
tion of the systems, and allows for a qualitative under- 
standing without the knowledge of the exact microscopic 
details. In this Letter, we analyze the appearance of cor- 
related manybody effects in a dense ensemble of atoms 
driven resonantly in a strongly interacting Rydberg state, 
and show that the behavior of the system can be under- 
stood in terms of a critical theory. 

Strongly interacting Rydberg atoms are an area of 
intense experimental investigation: the resonant diffu- 
sion of Rydberg excitations via the dipole-dipole inter- 
actions has been reported 0] , and the reduction of Ryd- 
berg excitation due to blockade effects has been observed 
d, 0, H, @, S Si- Furthermore, coherent optical exci- 
tation has been observed for individual Ryd berg atoms 
Q , essentially non- interacting ensembles 03, 1 1 lL , as 
well as in the strongly interacting regime [7| . Of special 
interest for the theoretical analysis are cold atomic sam- 
ples, where the thermal motion of the atoms is essentially 
frozen on the characteristic time scale of the Rydberg ex- 
citation ( "frozen Rydberg gas" ) Q ■ Then, the time evo- 
lution of such an atomic ensemble driven by a resonant 
Rydberg excitation has been extensively studied in the 
past using a numerical integration of the time-dependent 
Schrodinger equation for small sample sizes hj,hj and 
within a master equation approach [l5l . Hil ] . 

In contrast to the previous analyses studying the time 
evolution, the approach for the study of strongly inter- 
acting Rydberg gases presented in this Letter is based 
on the observation that the driven system relaxes into an 
equilibrium state. This equilibrium state is dominated by 
the thermodynamic phases of the Hamiltonian describing 
the driven system, which exhibits a continuous quantum 



phase transition in the detuning A; see Fig. 0Ja). Then, 
the resonant regime with A = is determined by the 
critical properties of the phase transition. We derive the 
universal exponents within mean-field theory, and show 
that the time evolution of the system and its relaxation 
into the equilibrium state is well described by a master 
equation derived from the microscopic Hamiltonian. 
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FIG. 1: (color online), (a) Phase diagram in the A-fi-plane: 
for a coupling Q — a second order phase transition appears 
from a crystalline phase (C) to a paramagnetic phase (PM), 
while at a detuning A = the system is dominated by the 
critical region, (b) Setup of the system with \g) the atomic 
ground state and |e) the excited Rydberg state coupled by 
driving lasers with Rabi frequency f2 and detuning A. 

We start with the Hamiltonian describing an ultracold 
gas of atoms driven into an excited Rydberg state with a 
repulsive van der Waals interaction. The relevant inter- 
nal structure for each atom is given by the atomic ground 
state \g)i and the excited Rydberg state |e)j, reducing the 
internal structure to a two-level system. The two inter- 
nal states are coherently coupled by external lasers with 
the Rabi frequency f2 and a detuning A; see Fig. [TJb). 
The characteristic time scale for a Rydberg excitation is 
short compared to the thermal motion of each atom, and 
the positions of the atoms are frozen : the positions 
rj are randomly distributed according to the distribution 
function of a thermal gas. After performing the rotat- 
ing wave approximation, the Hamiltonian in the rotating 
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frame can be written as a spin Hamiltonian (cf. 13]) 
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where P e e = l e )» ( e U = (1 + )/2 is the projector onto 
the excited Rydberg state of the i th atom. The last term 
accounts for the strong van der Waals repulsion between 
the Rydberg states with Cq oc n 11 , where n is the princi- 
ple quantum number of the Rydberg excitation [l3, EH . 
It is well known that dense samples of Rydberg atoms un- 
dergo collective ionization processes that might harm the 
coherent evolution [l^ . However for repulsive Rydberg S 
states these ionization processes are slowed down signif- 
icantly [20j ]. Recent experiments in very dense samples 
allowed to follow the coherent evolution involving Ryd- 
berg S states on timescales comparable to the timescales 
required to reach equilibrium • In these experiments 
the saturation is reached well before radiative or motional 
decoherence effects set in. Therefore the validity of the 
Hamiltonian ([I]) is well justified for an experimentally 
realistic situation. 

Within the standard experimental setup, the system is 
prepared with all the atoms in the ground state and the 
time evolution is studied by turning on the excitation 
lasers, i.e., \tp) = J\ i \g)i. While in absence of interac- 
tions, the time evolution of the Hamiltonian fl} results 
in coherent Rabi oscillations for each atom with the re- 
duced density matrix for each atom being a pure state, 
the interactions leads to correlations between different 
atoms, and eventually to a decoherence of this single 
atom pure state. Consequently, the density matrix de- 
scribing a subsystem of the atomic cloud will equilibrate 
into a steady state under time evolution as the surround- 
ing atomic states act as a reservoir interacting with the 
subsystem. This behavior is well confirmed by the exact 
numerical integration of the system with the Hamilto- 
nian fl}, see below, and is also observed in experiments 
via the saturation of the Rydberg excitation This 
observation opens up an alternative approach to study 
the long-time behavior of the system by focusing on the 
equilibrium states of the Hamiltonian fl}, and especially 
on its zero temperature phase diagram. 

For vanishing Rabi frequency il = 0, the exact ground 
state of the Hamiltonian JT]) can be analytically deter- 
mined in any dimension. It is dominated by a second 
order quantum phase transition for the critical detun- 
ing A c = 0; see Fig. [TJ For negative detuning A < 
the ground state is paramagnetic with all atoms in the 
atomic state \g)i) i.e., the experimentally relevant ini- 
tial state |V>) is the ground state of the system. In turn 
for positive detuning A > 0, the ground state prefers to 
excite atoms into the Rydberg state and the configura- 
tion minimizing the repulsive van der Waals interaction 
is obtained for a crystalline arrangement of the atoms. It 
is important to notice, that within the critical region the 



system is independent on the microscopic realization, and 
therefore, its properties are isotropic and homogeneous. 

The behavior of the system for finite f2 in the resonant 
regime A = is dominated by the critical behavior of the 
second order quantum phase transition. At the critical 
point with A = 0, the system is characterized by a single 
dimensionless parameter a = hfl/Cen 2 describing the 
ratio between the coupling energy Ml and the interaction 
energy Cqu 2 (here, n denotes the atomic density) . Then, 
the ground state properties such as the fraction of excited 
Rydberg atoms exhibit an algebraic behavior 



(2) 



with a universal scaling exponent v in the critical region 
with a < 1. In the following, we determine the critical 
exponent within mean-field theory. 

The important quantity in the mean-field theory is the 
averaged Rydberg fraction f R = (Pee ) . In addition, the 
van der Waals interaction gives rise to blockade phenom- 
ena; i.e., once a Rydberg state is excited the excitation of 
an additional Rydberg atom is strongly suppressed in the 
surrounding area characterized by a blockade radius a R . 
The correct description of this property is obtained by 
the pair-correlation function ^(r^ — Yj) = (Pel Pee )/f%, 
where 32 (r) vanishes in the blockaded region for |r| <C an, 
while at large distances |r| 3> a R the correlation dis- 
appears and (72 (r) = 1. The transition from a strong 
suppression to the uncorrelated regime is very sharp 
due to the van der Waals repulsion [13J, and the pair- 
correlation function is well described by a step function 
32(1") = Q (|r| — or). Then, the mean-field theory is ob- 
tained by replacing the microscopic interaction by the 
mean interaction of the surrounding atoms 

P£ P& « [P$ fn + fn ~ f%] 92 (r« - r,) , (3) 

which neglects the quadratic fluctuations around the 
mean field and reduces the Hamiltonian to a sum of single 
site Hamiltonians, H = J^-Hmf- In the scaling regime 
with a <C 1, the number of atoms in the blockaded regime 
is large, i.e., a R n 3> 1, which allows us to replace the 
summation over the surrounding atoms j by an integral 
over space with a homogeneous atomic density n. Then, 
we obtain the Hamiltonian for the i th atom, 
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T4p«4A=h-^+/,. (4) 
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with the characteristic energy scale Eq = Cqtl 2 . Note 
that in equilibrium the correlation function 172 (r) satis- 
fies the normalization condition nf R J dr [1 — 32(1")] = 1, 
which provides the relation a R — (S/infun) 1 / 3 . The ef- 
fective Hamiltonian is equivalent to a spin in a magnetic 
field h with h x = a/2 and h z = 87r 2 /^/9, and a constant 
energy offset ho = h z (l — f R ). Using a spin rotation, 
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FIG. 2: (color online). Numerical integration of the full 
Hamiltonian (crosses) for a — 1/36, N = 60 and averaged 
over 50 different initial conditions. The remaining oscilla- 
tions at large times are finite size effects. The solution of the 
master equation is shown as solid line. 
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we can diagonalize the Hamiltonian = ha z , 



with h = \Jh 2 x + h 2 z . Here, <j z % ) — cos 9a ( z ' + sin 9a 
denotes the Pauli matrix in the new basis with the rota- 
tion angle tan# = h x /h z . Within the new basis the sta- 
tionary equilibrium density matrix ffo) is diagonal. Its 
entries dependent on the experimental realization of the 
system: (i) for an adiabatic switching on of the Rabi 
frequency, the system remains in the ground state and 
the density matrix reduces to the lowest energy state 
of HmI- (ii) on the other hand, for a sudden switching 
on of the Rabi frequency, the density matrix is deter- 
mined by energy conservation: the energy of the equili- 
brated state is equal to the energy of the initial state, 

i.e., Tr|^F/5 (l) } = (ip\H\ip)/N = (here, N denotes 

the total number of particles). It can be checked that 
the scaling exponent remains the same in both cases, in- 
dicating that the equilibrium state in the second case 
is close to the ground state. In the following, we focus 
on a sudden switching on of the Rabi frequency. The 
mean-field solution for the equilibrium state reduces to 
(h /h)a z l ) /2, and a transformation into the 
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original coordinates yields the self-consistency relation 

(9a) 2 
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The solution in the limit a <C 1 provides the critical 
exponent v = 2/5 with the prefactor c = (9/16-7T 2 ) 2 / 5 . 

The same mean-field analysis can also be performed 
for arbitrary dimensions d < 4. The main modification 
is that the dimensionless parameter ad — hQ/Cen 6 ^ d ex- 
hibits a dependence on the dimension d and also h z obeys 
the modified scaling h z ~ /» ■ The general result pro- 
vides the scaling exponent — 2d/ (12 + d). 

We now compare the above mean-field solution to nu- 
merical studies of the time evolution with the full Hamil- 
tonian Eq. JTJ). We place N atoms randomly into a box of 
volume V having periodic boundary conditions. The di- 
mension of the full Hilbert space grows exponentially like 



2^; therefore, the exact dynamics can only be calculated 
for a relatively small number of atoms 13l.l2lj|. However, 
the strong van der Waals repulsion suppresses the occu- 
pation probabilities of many basis states, which allows us 
to significantly reduce the Hilbert space: for each basis 
state we compute the van der Waals energy and remove 
the state if its van der Waals energy is larger than a cutoff 
energy Ec- This reduction leads for N = 100 to approx- 
imately 10 6 relevant basis states compared to the 10 30 
basis states of the full Hilbert space. Convergence has 
been checked by increasing Ec- We have used a fourth- 
order Runge-Kutta scheme for performing the numeri- 
cal integration of the Schrodinger equation. The time- 
dependence of the Rydberg fraction /«(t) = Si(-^ee)/-^ 
is shown in Fig. [2] (crosses), for an average of 50 differ- 
ent uniformly distributed random initial conditions with 
N = 60 and a — 1/36. We find a clear saturation of the 
excited Rydberg fraction the remaining oscillations 
at large times are further suppressed for larger system 
sizes, and therefore represent a finite size effect. 

In order to investigate the scaling behavior of the satu- 
rated Rydberg fraction fn we have taken an average over 
the Rydberg fraction /.r(£) for times 250 < tC§n 2 jh < 
400 for 50 different random initial conditions. The di- 
mensionless parameter a has been varied by changing 
the number of atoms from N = 52 to N = 100, and 
by changing C§ from 0.01 to 0.04. The scaling behavior 
is shown in Fig. [3] The data show a power law depen- 
dence according to fa — ca v . The fit to the numerical 
data provides an exponent v — 0.404, which is in very 
good agreement with the result derived in the mean-field 
theory. For the onedimensional case the numerically ob- 
tained critical exponent is v 1D — 0.150. Surprisingly, this 
value is again very close to the mean-field prediction, 
indicating that the van der Waals interaction strongly 
suppresses quantum fluctuations. 



-1.0 



-2.0 



1 1 

3D 




v ss 0.404 






^6" 



-2.5 



-1.5 




-1.4 



-1.5 



-7.5 -7 

l°gl0 Q 1D 



FIG. 3: (color online). Numerical results for the saturated 
Rydberg fraction fn in a 3D and a ID setup: the system 
exhibits an algebraic behavior fa ~ a v with v ~ 0.404 and 
v 1D fa 0.150. 

Finally, we are interested in a description of the time 
evolution and derive a master equation with the mean- 
field solution as its stationary state. The natural mech- 
anisms for the equilibration into a stationary state are 
the residual interactions between the atoms, which go 
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beyond the mean-field description. We write the exact 
Hamiltonian {1} as a sum of the mean-field terms and the 
remaining fluctuations H = J2i + AH. The deriva- 
tion of the master equation from the microscopic Hamil- 
tonian uses the time-convolutionless projection operator 
method with an extended projection operator [22j |: we 
select a single site i, which will play the role of the sys- 
tem with the Hamiltonian H^l, while the surrounding 
atoms act as the bath coupled to the system state by 
the Hamiltonian AH. The role of the pair-correlation 
function is to enforce the blockade regime. Here, we 
assume the pair-correlation function to be fixed during 
the time evolution. Then, its influence is well accounted 
for by expressing the remaining interactions as AH = 



E»<j .92(r, - Tj) C 6 /\n - r,f (p« - f R ) (p« - f R 
The projection operator consistent with our mean-field 
theory, i.e., VAHp = reduces to Vp = i p^- Here, 
denotes the reduced density matrix defined by the 
partial trace p^ = Tr.; {p}, which performs the trace over 
all atomic states except of the i th atom. Next, it is useful 
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to express the operator P ee . , . 

of the projections onto the eigenstates of where 
fkuci = 2/i denotes the energy splitting of the mean-field 
Hamiltonian. Introducing the interaction picture with 



respect to J2i ^mf the interaction Hamiltonian H>ll (t) 
within the rotating wave approximation reduces to 

r(i,j) (,\ \ " ^6 



H, 



'(*)= E 
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The terms A± bJo describe the exchanges of an excitation 
between the system and the surrounding bath and is 
relevant for the equilibration of the time evolution, while 
the terms with Aq account for a dephasing. The second 
order time-convolutionless master equation 22, [23| for 
the reduces density matrix p^> (t) takes the form 
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with the rates 
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Here, we have again used the translation invariance of the 
system in the critical region a«l allowing us to replace 
the average over the j th atom, with the local density ma- 
trix p( l \t). The master equation is a highly nonlinear 
equation for the local density p^> (t) , where at each time 
step the mean- field f R = Tr{/>M pje^} and the rates 
have to be determined. We would like to stress that the 
master equation conserves energy with the above mean- 
field solution as a stationary state. The master equation 
can be efficiently solved numerically after a transforma- 
tion back into the Schrodinger picture. In Fig. [5] the so- 
lution of the master equation (solid line) is compared to 



the full dynamics (crosses). The pair-correlation function 
strongly influences the decoherence rates 7u),a/- 

As the 

pair correlation is not determined self-consistently within 
our approach, we apply once a single fit for the effec- 
tive Rabi frequency and the saturated Rydberg fraction 
to account for these corrections. Then, we find perfect 
agreement between the full dynamics and the solution of 
the master equation, see Fig. [2j the deviations at large 
times are accounted to the finite size effects of the full dy- 
namics. The characteristic time scale of the oscillations 
is determined by the collective Rabi frequency \/iV b i ock f2, 
with iV block ~ na R the number of atoms within the block- 
aded volume of the van der Waals interaction. 

In conclusion, we have established a universal scaling 
behavior in strongly interacting Rydberg gases due to the 
existence of a second order quantum phase transitions. It 
remains an open question, whether it is experimentally 
possible to reach the crystalline phase. 
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